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We consider a two-component Bose-Einstein condensate, which contains atoms with magnetic 
dipole moments aligned along the z direction (labeled as component 1) and nonmagnetic atoms 
(labeled as component 2). The problem is studied by means of exact numerical simulations. The 
effects of dipole-dipole interaction on phase separations are investigated. It is shown that, in the 
quasi-one-dimensional regime, the atoms in component 2 are squeezed out when the dimensionless 
dipolar strength parameter is small, whereas the atoms in component 1 are pushed out instead when 
the parameter is large. This is in contrast to the phenomena in the quasi-two-dimensional regime. 
These two components are each kicked out by the other in the quasi-one-dimensional regime and 
this phenomenon is discussed as well. 



I. INTRODUCTION 

Over the past decade, the experimental realization of a 
Bose-Einstein condensate (BEC) of 52 Cr atoms [IH1] with 
a strong dipole-dipole interaction has given new impetus 
to theoretical investigations of BECs with long-range in- 
teraction at low temperatures. Due to the large magnetic 
dipole moment of Cr atoms ncr = G^b (where \ib 

is the Bohr magneton), the anisotropic interaction be- 
tween polarized magnetic dipoles results in anisotropic 
deformations during expansion of the condensate. Dipo- 
lar quantum gases are governed by the d-wave symmetry 
of the long-range dipole-dipole interaction, which causes 
novel properties such as unusual stability properties [5], 
exotic ground states [5] [7], and modified excitation spec- 
tra 019]. For a review of dipolar condensates, see [ID] . 

The miscibility and separation of a two-component 
BEC have been studied theoretically [TTJ [T^] and ex- 
perimentally [13] previously. Recently, more detailed 
and controlled experimental results have been obtained, 
illustrating the effects of phase separation in a multi- 
component BEC [2]. In all these papers, the stud- 
ies of the binary condensates were limited to the case 
of s-wave interactions; while a great deal of attention 
has been drawn recently to dipolar BECs. Hexagonal, 
labyrinthine, solitonlike structures, and hysteretic behav- 
ior have been studied in a two-component dipolar BEC 
[15] . as well as immiscibility-miscibility transitions [16] . 

The theoretical description of a BEC in the dilute limit 
in the framework of the extended Gross-Pitaevskii equa- 
tion (GPE), which has been part of many text books on 
quantum mechanics or Bose-Einstein condensates [17] , is 
well known. The time-independent GPE reads 
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iP(r) 



(1) 



where \i is the chemical potential, m is the atomic 
mass, n(r) = iV|-0(r)| 2 is the density of atoms, 
which is normalized via J n(r)d 3 r = N, and V(r) = 
m (oj 2 x 2 + w 2 y 2 + uj 2 z 2 ) /2 is the external harmonic trap 
confining the gas. The local mean-field potential 
U\i/j(r)\ 2 represents the s-wave interaction, with U = 
4nh 2 a s /m and scattering length a s . The long-range in- 
teraction between two particles, $ dd (r), represents the 
dipolar interaction, which is given by [3 [8] 



$dd(r) = Jd 3 r'U dd (r-r')\yj{r') 



(2) 



Here U d d{r) = [C dd / (An)} e.^j (<% - 3fify) /r 3 describes 
the interaction between two dipoles which are aligned by 
an external field along a unit vector e through a distance 
r. The coupling C dd = E 2 a 2 /en where the dipoles are 
induced by an electronic field E = Ee, with the static 
polarizability a and the permittivity of free space eo- If 
the atoms have a magnetic dipole moment d m aligned by 
a magnetic field B = Be, the coupling is C dd = Hod^ 
with the permeability of free space /Zo [IS]. A measure- 
ment of the strength of the dipole-dipole interaction rel- 
ative to the s-wave scattering energy is provided by a di- 
mensionless quantity (the so-called dimensionless dipole- 
dipole strength parameter) e dd = C dd / (3£7). For dipoles 
aligned along z, this dipolar mean- field potential ([2| can 
be expressed in terms of a fictitious electrostatic potential 
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where n(r) and <j>{r) satisfy Poisson's equation V 2 = 
—n(r). This formulation of the problem allows us to im- 
mediately identify some generic features of dipolar gases. 
For example, if n(r) and hence 4>(r) is uniform along the 
polarization direction z, then the nonlocal part of the 
dipolar interaction vanishes because of the operator d 2 . 
The details are in Sec lIIII 

The two-component dipolar BEC, confined in a cylin- 
drical trap, is described by two coupled Gross-Pitaevskii 
equations. We take Cr as component 1 and Rb as com- 
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ponent 2, then the GP equations can be written as 
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where Vi(r) 



/2 (i = 1, 2) represents 



a cylindrical harmonic trap with radial trap frequency 
LUip and axial trap frequency iOi Z for the ith component, 

and p = \J x 1 + y 2 . ij)i(r) = 4'i{p)4 , i( z ) ls the wave func- 
tion of component i with particle number iV,-. m; and 
lii are the mass and the chemical potential for the ith 
component, respectively. 

The interatomic and the intercomponent s-wave scat- 
tering interactions are described by Ui and U\ 2 , respec- 
tively, with the following expressions |17j : 
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where is the scattering length of component i and a\ 2 
is that between component 1 and 2. 

The goal of this paper is to analyze the properties 
of phase separation in a two-component dipolar BEC 
in quasi-one and quasi- two dimensions, and it is or- 
ganized as follows. In Sec. [IlJ the Gross-Pitaevskii 
equations of the two-component dipolar BEC are re- 
duced to quasi-one-dimensional form. Analytical insights 
into the phenomenology of phase separation in the ax- 
ial direction are shown. In Sec. Ill the quasi-two- 
dimensional Gross-Pitaevskii equations representing the 
two-component dipolar BEC are derived, and analytical 
results for the phase separation in the radial direction are 



given. Finally, the paper is concluded in Sec IV 



and (|6| as two dimensionless coupled differential equa- 
tions: 
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where a m — m 2 /mi and a w = lu 2z /uji z . 
Using the finite-difference approximation 
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we can write Hi as a symmetric tridiagonal matrix with 
diagonal elements v\ and subdiagonal or superdiagonal 
elements w\ in Hi, where 
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with mesh length h in the z direction. Here, free bound- 
ary conditions should be applied as lim ijji(z) = 0. Di- 

2— f ±00 

agonalizing these two symmetric tridiagonal matrices, we 
can obtain the ground state wave functions tpi (z) and the 
density profiles of the dipolar BEC [T2] • 

In this case, the critical value af 2 which causes phase 
separation is 



II. QUASI-ONE-DIMENSIONAL REGIME 
(INFINITE PANCAKE) 

We begin by neglecting the radial trapping frequency 
(uiip = 0) and increase the axial trapping frequency (uj iz ) 
to the extent that a highly flattened pancake-shaped 
BEC is produced with infinite radial extent. In this 
case, the trap has infinite radial extent and uniform two- 
dimensional (2D) density in the p direction. The Poisson 



equation V 



-n(r) reduces to d 2 



([3]) reduces to the contactlike form [TU] 



(2;), and Eq. 
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If we let yjfi/ (toiWi z ) and Tiu>i z be the units for length 
and energy, respectively, then we can rewrite Eqs. ^ 
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arising from U~l U 2 — U 2 2 = with Ul — 
Ui (1 + 2e dd ). This criterion for the strength of the re- 
pulsive intercomponent interaction is independent of the 
atom numbers as well as of the trap strength [21] . 

In our calculations, we take the scattering lengths of 
52 Cr and 87 Rb as 5 nm [3] and 10 nm [22], respectively. 
For the trap, we assume uj\ z = 2irx 160 Hz and a m a 2 = 1. 
In this case, we take the particle number as N% = N 2 = 
1000. 

First, we assume e dd = —0.3 and calculate the density 
profiles of those two components for various a\ 2 using 
the coupled GP equations. For a relatively small value of 
the intercomponent scattering length a 12 (4 nm in Table 
|l|), the two-component BEC does not separate [see Fig. 
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TABLE I. Dimensionless values of interactions in quasi-one 
dimension 
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[l] (A)]. Then we increase a%2 to the critical point (4.33 
nm in Table [I]), the phase separation of Cr and Rb atoms 
occurs and the density of Cr increases [see Fig. [I] (B)]. 
The dipoles in this case lie predominantly side by side, 
and the net dipolar interaction is attractive (repulsive) 
when Sdd < (sad > 0) [19]. As a result, the effective 
interatomic interaction U^^ 1 (containing the dipolar in- 
teraction and s-wave interaction) in Cr is less than the 
interatomic interaction U2 in Rb, and both of them as 
well as the intercomponent interaction are repulsive (see 
Table |l| . Then Cr pushes Rb out toward the edges of 
the trap, while the rubidium atoms also "squeeze" the 
chromium atoms, i.e., they act like a trap which enhances 
the trapping force on the chromium atoms and makes the 
peak density of Cr higher. In the one-dimensional case 
containing only s-wave scattering interactions, the crite- 
rion Eq. (10) works perfectly [23]. However, when we 
increase eJcTT as in Fig. [5] where e d d — 0.1), the antici- 
pated phase separation doesn't occur at the critical point. 
Instead, the two components admix homogeneously [see 
Fig. [2] (A)]. By calculating the interactions' strengths, we 
know that the two interatomic interactions and the in- 
tercomponent interaction have the same value (see Table 
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FIG. 2. (Color online) Ground-state density profiles in the 
two-component dipolar BEC for Edd = 0.1 and different values 
of the intercomponent scattering length: (A) ai2 = 7.5 nm 
(the critical point); (B) 012 = 7.6 nm. 



[T]); thus the two-component Bose gas behaves like a one- 
component Bose gas since all particles have, so to speak, 
the same hard-core diameter. As a result, all the atoms 
are in a homogeneous mixture. Increasing a±2 slightly 
beyond the critical value of the intercomponent scatter- 
ing length, we find that the two components are repelled 
by each other. Since the intercomponent interaction is 
larger than the two interatomic interactions as shown in 
Table |l] none of the atoms stays at the center of the trap 
[see Fig. |j(B)]. 

Increasing Edd (as in the Fig. [3] where Edd — 0.5), 
we find that the phase separation is different from that 
shown in Fig. [T] At the critical point, the rubidium 
atoms still occupy the center of the trap and their peak 
density becomes higher, while the chromium atoms are 
pushed out instead [see Fig. [3] (B)]. As we can see, 
this is because the effective interatomic interaction of the 



Q 






z(nm) 



FIG. 



1. (Color online) Ground-state density profiles in the 
two-component dipolar BEC for Edd = —0.3 and different val- 
ues of intercomponent scattering length: (A) 0,12 = 4 nm; (B) 
012 = 4.33 nm (the critical point). In this figure and Figs. [2] 
and [3] the blue solid line represents the density profile of 52 Cr 
and the red dashed line is that of 87 Rb. 
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FIG. 3. (Color online) Ground-state density profiles in the 
two-component dipolar BEC for Edd = 0.5 and different values 
of the intercomponent scattering length: (A) ai2 = 8 nm; (B) 
ai2 = 9.68 nm (the critical point). 
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chromium BEC is now stronger than the interaction of 
the rubidium BEC (see Table [I]), and the rubidium atoms 
get an extra trapping force from the chromium atoms. 
Based on this, we may control the component which we 
intend to move out of the trap center to some extent. 



III. QUASI- TWO-DIMENSIONAL REGIME 
(INFINITE CIGAR) 



We now consider the two-component dipolar cigar- 
shaped BEC in a quasi-two-dimensional regime using a 
similar methodology to that for the pancake shape. If 
we neglect the axial trapping (u>^ z = 0) and consider the 
BEC to be uniform along z, then Eq. ^ reduces to the 
contactlike form 1191: 



$dd = -U x e dd n (p) = -E/iJVieddl^i (p) | 2 . (11) 
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with mesh length h in the p direction. By diagonalizing 
these two nonsymmetric tridiagonal matrices, we can get 
the ground state wave functions ^ (p) and the density 
profiles of the dipolar BEC [IS]. Here, the critical value 
a\ 2 is obtained by 



Then we can rewrite Eq. ^ and |6]) as two dimensionless 
coupled differential equations: 
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We assume that u>i p = 2irx 160 Hz, the particle number 
is taken as Ni = No. = 2000, and d ro o„ = 1 as well. 
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with units of length and energy and a m = m 2 jm\ as in 
Sec. |Hj while a w = ui 2p /oji p is just different in form from 
that in Sec. [TTJ Here, for simplicity and comparability, we 
assume that the values of these two a w are equivalent, i.e., 
w 2z /wi 2 = uj 2p /u lp . 

Using a similar methodology, we apply boundary con- 
ditions lim ipi(p) = and -i~ipi(p)\ P =o — 0, and the 
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finite-difference approximations |20j 
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then we can write H t as a nonsymmetric tridiagonal 
matrix containing diagonal elements sj, subdiagonal ele- 
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FIG. 4. (Color online) Ground state density profiles in the 
two-component dipolar BEC for Edd = —0.3. The intercom- 
ponent scattering length a\2 = 5.5 nm in (A) and (B), while 
«i2 = 7.81 nm (the critical point) in (C) and (D). In this fig- 
ure and Figs. [5] and [6] (A) and (C) represent the chromium 
atoms, while (B) and (D) represent the rubidium atoms. The 
numerical values of the units for density and length are re- 
spectively 1 x 10 8 cm -2 and 1 /im in our calculations. The 
field of view for (A)-(D) is 10 x 10 ^m. 



TABLE II. Dimensionless values of interactions in quasi-two 
dimension 
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Jjeff2 (1Q - 3) ^ (1Q - 3) jj 2 (1Q -3) 



-0.3 
0.1 
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82 



57 



31 



55 (ai2 = 5.5 nm) 
78 (af 2 = 7.81 nm) 
30 (ai2 = 3 nm) 
65 (012 = 6.5 nm) 
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49 (afg = 4.84 nm) 
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In this case, the dipoles are predominantly end to end, 
and the net dipolar interaction is repulsive (attractive) 
when Edd < (Edd > 0) [15] . This is in contrast to 
the case in the quasi-one-dimcnsional regime. We then 
find that the chromium atoms are squeezed out when 
£dd = —0.3 at the critical point, and form a shell around 
the rubidium atoms which are situated at the center of 
the trap [see Fig. [4] (C) and (D)]. The critical point 
of di2 in this case is not the same as in the quasi-one- 
dimensional regime, but higher. We believe the cause 
is that the effective interatomic interaction {/^ * 2 of Cr, 
which is repulsive, becomes stronger than in the former 
case (see Tables [I] and . We cannot observe a similar 
phase-separated state for Edd = 0.1 as in Fig. [2](B), other 
than in the quasi-one-dimcnsional regime. As shown in 
Fig. [5] the two components are mixed homogeneously 
at a small value of a\2 (3 nm in Table |n]), and phase 
separation occurs at the critical point. 

On increasing Edd to 0.5, we find that the attractive 
strength of the dipolar interaction becomes very large, 
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FIG. 6. (Color online) Ground state density profiles in the 
two-component dipolar BEC for Edd = 0.5. The intercom- 
ponent scattering length a± 2 = 2 nm in (A) and (B), while 
ai2 = 4.84 nm, the critical point, in (C) and (D). The units 
are the same as in Fig. [4] The field of view for (A)-(D) is 
10 x 10 fim. 



and the chromium atoms begin to occupy the center of 
the trap, while the rubidium atoms form a shell around 
chromium (see Fig. |6j. The phase-separated states for 
the same Edd are totally different between the quasi-one- 
and quasi-two-dimensional regimes. The interactions are 
known to be the predominant factors which can affect the 
ground-state density profile of a BEC, so different prop- 
erties of the dipolar interaction can make a difference to 
the phase separations between the two regimes. For the 
same value of Edd, the repulsion or attraction of the net 
dipolar interaction is opposite in the two regimes, so the 
effective interatomic interaction of Cr is totally differ- 
ent between the two regimes (compare Uf f2 of Table [□] 
with Ul^ 1 of Table H. If Edd is negative, the repulsion 
of the effective interatomic interaction in the quasi-two- 
dimensional regime is much stronger than in the quasi- 
one-dimensional regime (see Tables [I] and [TT]) , so the ru- 
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FIG. 5. (Color online) Ground state density profiles in the 
two-component dipolar BEC for Edd = 0.1. The intercom- 
ponent scattering length a\ 2 = 3 nm in (A) and (B), while 
ai2 = 6.5 nm, the critical point, in (C) and (D). The units 
are the same as in Fig. [4] The field of view for (A)-(D) is 
10 x 10 Aim. 
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FIG. 7. (Color online) Ground state density profiles of com- 
ponent 1 (Cr). Ni = 2 x 10 3 , Edd = -0.3, a 12 = a\ 2 = 7.81. 
From (A) to (C), N 2 = 2 x 10 3 , N 2 = 2 x 10 4 , iV 2 = 2 x 10 5 , 
respectively. The units are the same as in Fig. [4] The field 
of view for (A)-(C) is 20 x 20 fim. 
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FIG. 8. (Color online) Ground state density profiles of com- 
ponent 2 (Rb). N 2 = 2 x 10 3 , e dd = -0.3, aia = a\ 2 = 7.81. 
From (A) to (C), Ni = 2 x 10 3 , Ni = 2 x 10 4 , Ni = 2 x 10 s , 
respectively. The units are the same as in Fig. [4] The field 
of view for (A)-(C) is 20 x 20 /im. 

bidium atoms are held at the center of the trap while 
the chromium atoms are forced to form a shell around 
the center at the critical point of phase separation. The 
situation for a relatively large positive Edd is the opposite. 

In addition to the strength of the interactions, the den- 
sity profiles also depend on the particle numbers iVi and 
7V2 [12] . Thus we do some calculations where we hold con- 
stant the number of one component's atoms and change 
the other in the quasi-two-dimensional regime. Clearly, 
the number of atoms plays a key role in the phase separa- 
tion of a two-component BEC. In the case of Edd = —0.3, 
Cr forms a shell around Rb. When the number of Cr 
atoms is fixed, they are squeezed further out as the num- 
ber of Rb atoms increases (See Fig. [7]). If we fix the 
number of Rb atoms instead, one can see that the rubid- 
ium atoms are compressed further as the number of Cr 
increases (See Fig. [8]). Finally, we would like to mention 
that the situation of the quasi-one-dimensional regime is 
similar to this case herein. 



IV. CONCLUSION 

In the present work, we have analyzed the phase sep- 
aration of a two-component dipolar BEC in the quasi- 
one- and quasi-two-dimensional regimes, as a natural pa- 



rameter of the system (the intercomponent interaction 
strength) was varied. Our analysis was presented for 
the case of a two-component Bose-Einstein condensate 
trapped in a harmonic-oscillator potential. 

We reduced the long-range dipole-dipole interaction to 
a contactlike form, when a dipole moment is aligned along 
the z direction, and then used a numerical method to 
solve two coupled nonlinear Gross-Pitaevskii equations in 
the low-dimensional regime. We were able to elucidate 
the phase separation condition for the intercomponent 
coupling parameter as the dimensionless dipolar interac- 
tion strength parameter Edd was varied. 

Both repulsive and attractive dipolar interactions have 
been taken into consideration. Among all the phase- 
separated states, Edd plays a significant role. By mod- 
ulation of its strength, the two components might alter- 
nately stay in or be forced to be outside the center of 
the trap. In addition to the interactions, we considered 
the dependence of the density profiles on the number of 
atoms as well. 

We also described the differences between the quasi- 
one- and quasi-two-dimensional regimes in detail. Due to 
the anisotropic property of the dipole-dipole interaction, 
phase separation in the axial direction differs from that 
in the radial direction. 

A natural extension of this work is to try to gener- 
alize the ideas presented herein to higher dimensions or 
larger number of components (i.e., spinors). Especially in 
the former case, more complicated pattern formations on 
the interface, such as Rosensweig hexagonal peaks and 
a labyrinthine pattern, have been observed [15] . Such 
studies are presently in progress and will be presented in 
future work. 
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